LU TP 99-07 
June 1999 



On Particle Multiplicities in Three-jet Events 



Patrik EdenQ, 
Gosta Gustafson-'^, 
Valery KhozeQ 



Abstract: 

A thorough verification of the distinct differences in the properties of quark and gluon jets 
is considered as one of the most instructive tests of the basic ideas of QCD. In the real life 
experiments such a comparison appears to be quite a dehcate task and various subtle issues 
require further theoretical efforts. In this paper we discuss in detail the possibility to extract 
the theoretically adequate information from the particle multiplicity patterns in three-jet 
events in e'^e~ annihilation. 
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1 Introduction 



As well known, the larger colour charge of gluons (C^ = A^c = 3) compared to quarks 
(Cp = (Ar,2 _ l)/2iVc = 4/3) leads to various distinctive differences between the two types 
of jets, for recent articles see e.g. and the review |^. Thus, a detailed comparison of the 
properties of quark and gluon jets provides one of the most instructive tests of the basic ideas 
of QCD. An experimental verification of these differences has been a subject of quite intensive 
investigations, especially in the last years, e.g. 0. However, obtaining of the theoretically 
adequate information about the properties of the gluon jet appears to be not an easy task. 
Recall that the analytical QCD results address the comparison between the energetic gluon 
and quark jets emerging from the point-like colourless sources, and that (unlike the qq case) 
the pure high energy gg events at present are not available experimentally.^ 

So far, most studies of the structure of gluon jets have been performed in three-jet events of 
e'^e~ annihilation. As a rule, these rely on a jet finding procedure both for selection of the 
qqg events and for a separation between the jets in an event. Without special care, such an 
analysis is inherently ambiguous and may suffer from the lack of the direct correspondence 
to the underlying theory. Recently some more sophisticated approaches have been exploited 
(see e.g. [3, 5-8]) which allow better theoretical significance. There are still a number of 
issues which are frequently overlooked in the present gluon jet analyses and some further 
theoretical efforts are required. First of all, this concerns particle multiplicity distributions in 
the jets. Clarification of these issues is the main aim of this paper. More detailed description 
of the theoretical framework can be found in ref. [0. 

In particular the following problems are addressed. 

1. Different approaches to the three-jet studies employ different definitions of the qqg kine- 
matics. In particular, this concerns such a key variable as a transverse momentum scale of 
the gluon, p±. Our first issue here is to discuss an exact definition of this quantity, which 
governs radiation from the gluon. 

2. The definition of the three-jet topology with the gluon registered at a given p± imposes 
an obvious requirement that there are no other subjets in the event with the transverse mo- 
mentum exceeding p±. We have to investigate quantitatively the impact of this requirement 
on the jet sample. 

3. To calculate predictions from perturbative QCD, using the assumption of local parton 
hadron duality (LPHD) 0, a cutoff is needed for the infrared singularities. As discussed 
in detail in ref. |I| such a cutoff depends on the soft hadronization process and can not 
be uniquely specified from perturbative QCD alone. Thus, the result is necessarily model 
dependent. 

In what follows we discuss these three issues successively in sections 2, 3 and 4, and in section 
5 we study their effect on analyses of 3-jet events in e"'"e~-annihilation. 

•^In principle, it is possible to create a pure source of the colour singlet gg events at a future linear 
e+e~ collider through the process 77 ^ gg ||]. 
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2 Definition of p_i 



In the simplest case of soft radiation, p± can be easily defined, as the quark and antiquark 
specify a unique direction. For large p± gluons, however, the q and q get recoils such that 
there is no obvious direction against which the transverse momentum should be measured. 
To have well defined expressions such a direction has to be specified. In the Lund dipole 
formalism [10-16] p± has been defined according to (subscript Lu for Lund) 

pIlu - (1) 

where Sqg denotes the squared mass of the quark-gluon system etc. In this particular frame 
the gluon rapidity is given by the expression 

l/ = ^ln(^). (2) 
The kinematically allowed region is given by 

P±Lu < \y\ < In (^^) = ^(L - «:lu); L ^ In(^), n^^ ^ H^)- (3) 

These variables have the advantage that the phase space element usually expressed in the 
scaled energy variables Xq and Xq is exactly given by the simple relation 

sdxqdxq = dp^Ludy- (4) 
As discussed in section ^, p_llu may also work well as a scale parameter in the QCD cascade. 



An alternative definition has also been used in the literature, e.g. by the Leningrad group [17 
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piLe - (5) 
■5qq 

This definition corresponds to the gluon transverse momentum in the qq cms (with respect 
to the qq direction). It is notable that in this frame the gluon rapidity is also exactly given 
by the expression in Eq @. The two p^- definitions agree for soft gluons, but deviate for 
harder gluons. While p^Lu is always bounded by y/s/2, p±he has no kinematic upper limit 
in the massless case. 



3 Bias from restrictions on subjet transverse momenta 



The effect of a cutoff in p± has been discussed previously [[T^, |19| . Here we give a brief review 
of the results, in order to end the section with an investigation of the numerical importance 
of subleading terms. These are essential for a correct analysis of three-jet events, which will 
be discussed in section |[ 
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To see the qualitative features of the bias we first study e~^e~ qq events within the Leading 
Log approximation (LLA). The quark and antiquark emit gluons according to the well-known 
radiation pattern 



dng ^ Cy—-. -7-^ = Cf —<^y = Cp My; k = ln(p^/A ). (6) 

TT (1 - Xq)(l - ,Tq) 7^ P± 71" 



We have here used Eq (Q), and in the following we define p± and y according to Eqs ([l|) 
and (0), unless otherwise stated. 

Due to colour coherence the hadronic multiplicity Ng{K) in a gluon jet depends on the p± 
of the gluon and not on its energy (see, e.g., refs [HI)- Summing up the contributions 
from all gluons in a cascade we arrive at the average multiplicity iVqq(L = ln(s/A^)) in the 
original qq system [13, 15-18] (Refs [15-18] include also nonleading terms.) 



A^qq(L) ^ Tdfi: /'^"^ dyCp^^N^iK) = ("^ <1k{L - t^)C^'^^N^{K). (7) 



(We have here introduced a lower cutoff kq for the integral over transverse momentum. This 
point will be discussed in section ^.) Taking the derivative with respect to L we find 
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Consider now a sample of events selected in such a way that there are no subjets with 
P± > P±cut- (Within a k±-hased cluster scheme with a resolution parameter p_Lcut, this 
means that there are only two primary q and q jets.) To obtain the multiplicity A'"qq(L, Hcut) 
in this biased sample, we must restrict the k integral in Eq (^) to the region k < Kent- We 



then find [11 



Nq^{L, Kcut) ~ A/'qq(/tcut) + (L - K^nt) %qif^cnt) ■ (9) 

The first term corresponds to two cones around the q and q jet directions. Here the p± 
of the emissions is limited by the kinematical constraint in Eq (^ rather than by Kent- It 
also corresponds exactly to an unbiased qq system with cms energy p±cut- The second term 
describes a central rapidity plateau of width (L — K^nt), in which the limit for gluon emission 
is given by the constraint Kent- This expression for a two-jet event can be generalized for 
a biased multi-jet configuration, and a similar discussion applies also to the multiplicity 
variance, cf. ref |TT[. (Similar equations for biased two-jet and three-jet events were later 



discussed also in ref |T9|] .) 

The average particle multiplicity in the selected two-jet sample is smaller than in an unbiased 
sample. The modification due to the bias is similar to the suppression from a Sudakov form 
factor. It is formally 0{as), but it also contains a factor ln^(s/p^). Thus, it is small for 
large p_L-values but it becomes significant for smaller p±. This clearly demonstrates that 
the multiplicity in this restricted case depends on two scales, a/s and p±cut- The p± of an 
emitted gluon is related to the virtual mass of the radiating parent quark. Therefore, the 
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Figure 1: The effect from the bias due to a constraint p±cut on emitted suhjets, at 90GeV 
energy. The figure shows the ratio of biased over unbiased multiplicities as a function ofp±cut ■ 
The results for LLA and MLLA relations (Eqs (j^ and (jl^, respectively) differ significantly 
from each other. The result of the MLLA relation in Eq pZ[j, using the p± definition in 
Eg is in good agreement with Ariadne MC and Durham cluster algorithm results. 



two scales ^/s/2 and p±cnt represent the energy and virtuality of the quark and antiquark 
initiating the jets. 

Though the LLA result in Eq (|) describes the qualitative features of the bias, subleading 
corrections are needed for a quantitative analysis. Within the Modified Leading Log ap- 
proximation (MLLA) ||20[ , subleading terms are included, which affect the prediction for the 
unbiased multiplicities and, thus, implicitly also the biased multiplicity in Eq @. Further- 
more, it is in [0] shown that the expression in Eq for the biased multiplicity is explicitly 
changed when MLLA corrections are considered. An unbiased system should be restored 
when p±cnt approaches the kinematical limit ^/i/2, but the r.h.s. of Eq (H) equals the un- 
biased quantity Aqq(L) only when p±cut = y/s. The relation consistent with the MLLA 
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Aqq(L, Kcut) ~ Aqq(Kcut + Cq) + (L - Kcut " Cq) Aq-(Kcut + Cq) ; Cq = - . (10) 



The bias is illustrated in Fig The dotted line shows results from the Ariadne MC [0] , 
when the Durham cluster algorithm [^] is used to define a biased sample of events classified 



as two-jet events with a ?/cut equal to p±cut/^- The MC results agree well with the prediction 
of Eq (p!OD, where for p±cut we have used the p_L-definition in Eq (|l|) (solid line). The predicted 
effect is below 5% for pxcut > 20GeV, but increases rapidly for smaller p_Lcut- 

Fig m presents also the result using the LLA relation in Eq (|) (dashed hne). To elucidate 
the effect of the differences between Eq (^ and (^^, we have used the same expression 
for the unbiased quantities Aqq and N^^. (These are obtained by a simple fit to Ariadne 
MC results, which are in good agreement with the MLLA.) As seen, the subleading terms 
are important; the LLA relation significantly overestimates the effect. To our knowledge 
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experimental data for this bias have not been presented. Such data should be obtainable 
in a rather straightforward analysis, which thus readily could test the accuracy of the MC 
result or the MLLA relation. 



4 Infrared cutoffs 



Gluon radiation diverges for coUinear and soft emissions. Therefore, to estimate the hadronic 
multiplicity from the assumption of LPHD |^, a cutoff is needed. Naturally, the cutoff must 
be Lorentz invariant. For collinear emissions a single Feynman diagram dominates, and 
there are two possibilities, the virtual mass, /i, of the emitting parent parton or the transverse 
momentum, p^, of the emitted gluon measured relative to the parent parton direction. These 
quantities are connected by the relation 



pi = ^^'z{l-z), (11) 

where z equals the light cone momentum fraction taken by the emitted gluon. The transverse 
momentum is directly related to the formation time, and, therefore, we regard this as the 
most natural choice for a cutoff. (For a further discussion see ref |T[|.) 

For soft emissions no obvious cutoff is available, however. As several Feynman diagrams 
contribute and interfere, there is no unique parent parton. Consequently /i^ or cannot be 
uniquely specified and, therefore, cannot be directly used. (Obviously a cut in energy is not 
possible, as this is not Lorentz invariant.) 

For soft emissions from a single qq colour dipole a cutoff in p_i_ is still the natural choice 
if measured in the cms, where the q and q move back to back. For emissions from a more 
complicated state the situation simplifies greatly in the large- A^c limit, as many interference 
terms disappear. In this limit the emission corresponds to a set of independent colour 
dipoles p2|, |12|. The natural choice for the cutoff is then p± in the cms of the emitting 
dipole (measured with respect to the dipole direction). We note that this implies that the 
soft gluons connect the hard partons in exactly the same way as the string in the string 
fragmentation model ||23|, which illustrates the connection between perturbative QCD and 
the string model [1 



For the physical case with 3 colours, extra interference terms appear with relative magnitude 
1/A^ O) HH • Here nonplanar Feynman diagrams contribute, and it is impossible to uniquely 
specify a parent parton or a relevant p±. Thus, a more fundamental understanding of 
confinement is needed to specify the cutoff, which cannot be determined from perturbative 
QCD alone |l|. In hadronization models the interference terms correspond to the 

problem of "colour reconnection" , and different models have been proposed ||2^. None of 
these can be motivated from first principles, and only experimental data can differentiate 
among the various models. 



In spite of the formal uncertainties, the success of current Monte Carlo programs ^ 
indicate that the colour suppressed interference terms do not have a very large effect. This 
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is also supported by recent searches by OPAL of the reconnection effects in hadronic Z 
events 0]. In most parton cascade formahsms, a cascade cutoff motivated in the large- A^c 
limit is used also for finite A*",,. The colour interference effects are accounted for by reducing 
the colour factor from A'c/2 to Cp in regions coUinear with quarks and antiquarks, and, due 
to colour coherence, also in some parts of the central rapidity region. We note, however, 
that some subtle interference phenomena, as a matter of principle, cannot be absorbed into 
a probabilistic scheme, see for details. These are still awaiting a thorough experimental 
test. 



5 Formalism for three-jet events 



After these general discussions we are now ready to consider three-jet qqg systems. To 
simplify the discussion we first study the large- A'c limit. The emission of softer gluons from 
a qqg system corresponds then to two dipoles which emit gluons independently. If a gluon 
jet is resolved with transverse momentum this imposes a constraint on the emission of 
subjets from the two dipoles. Thus, the contribution from each dipole is determined by an 
expression like Eq (p^Ol). For relatively soft primary gluons the constraint should be given 
by p±cut = P±g- For hard gluons p_llu is of the same order as its parent quark virtuality, 
and in ref it is shown that C(a;^) matrix elements are well described if p_llu is used as 
an ordering parameter for the perturbative cascade. This is also indicated by the successful 
applications of the Ariadne MC. We will, therefore, assume that the constraint on further 
emissions is well described by the identification p_Lcut = P±Lu- The multiplicity in a qg dipole 
with an upper limit on -p^ can, just as for the qq case discussed in section |^, be described as 
two forward jet regions and a central plateau. 

We note that if the three-jet events were selected using a cluster algorithm with a fixed 
resolution scale , then the constraint on subjet transverse momenta, p^cut? would be smaller 
than the of the gluon jet (as the gluon jet was resolved). In this case most jet definitions 



give three jets which are all biased ||TT|, |19[. We will, however, here focus on three-jet 



configurations obtained by iterative clustering until exactly three jets remain, without a 
specified resolution scale, where hence the constraint on subjet p± is described by pxcut = 
Pllu- As we will see, this implies that the bias on the gluon jet is negligible, which makes 
this selection procedure suitable for an investigation of unbiased gluon jets. 

For finite Ac the different dipoles in a multi-parton configuration can not be completely 
independent of each other. However, encouraged by the success of MC programs, let us 
assume that the main effect of finite A'c is that the colour factor, which determines softer 
gluon emission, is reduced from Nc/2 to Cp in the domains where the emission is dominated 
by radiation from the quark or the antiquark leg. Let us assume that a rapidity range Yq in 
the qg dipole is similar to a corresponding range in a qq dipole, while the remaining range 
Lqg — Yq is similar to a range in one half of a gg system. The corresponding ranges in the 
gq dipole are Yw and Lgq — Y^. This implies that the total multiplicity in the qqg event 
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corresponds to the expression 



For the constraint p^cut we have here written klu, which is appropriate for the selection 
procedure discussed above. 

As discussed in section 4, the size of Yq and Yq cannot be uniquely determined within 
perturbative QCD. Possibly the most natural choice is to assume that the quantity Yq + Y^ 
corresponds to the energy in the qq subsystem [|18|, which implies 



Yq + Yq^ ln(Sqq/A2) = Lqq. (I3a) 

The relation in Eq (|T3|a) can be regarded as an educated guess, but a finite shift cannot be 



excluded. In ref [jT5[ it is assumed that 

Yq + Yq^ Hs/A^) = L, (13b) 



which agrees with Eq ([13| a) to leading order. For relatively soft gluons we have Sqq ^ s, and 
in this case Eqs (|13|a) and (|13|b) are approximately equivalent. The assumption in Eq (0a) 
implies that the energy scale for the gluon term is given by Lqg + Lgq — Lqq = klc- Similarly 
we get from Eq (|13|b) the corresponding gluonic energy scale k^^. 

The effect of the p± constraint is rather different in the two terms in Eq (|T2p. For the 
gluon term the energy scale is in general only slightly larger than the bias scale klu- This 
implies that in most cases the bias can be disregarded in this term. Inserting the different 
assumptions in Eqs ([T3|a) and (0b) into Eq (|12]) then gives 

A'qqg ^ Nqq{Lqq, Km) + ^^"§§(^10), (14a) 
iVqqg ~ Nqq{L, ^Lu) + ^N^git^l^u) ■ (14b) 

We note that the consistency between Eqs (p!^) and (0b) follows from the fact that the 
total rapidity range in the two dipoles, Lqg + Lgq, can be expressed in two different ways by 
the equalities Lqg + Lgq = Lqq + = L + klu. In particular, we see from these equalities 
that the argument in iVgg has to be p^Le ^'i dH^) ^^^1 Pj^lu -^1 (p^b), and not e.g. 

The leading effect of a finite shift inYq + Yq is colour-suppressed, and therefore not expected 
to be large. However, subleading corrections introduce a difference between the results of 
Eqs ([T%i) and (^4|d). This is seen in Fig ^, where the difference is approximately 1 particle 
for y/s^ = 60GeV. In the calculations of A'qqg in Fig we have used the expressions in [|T| 
for the multiplicities A'^qq and Agg. These include MLLA corrections and recoil effects, which 
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Figure 2: Nq^g as a function of p^i,^ for 
y/s^ = 60 GeV. The different predictions 
from Eqs (j7^a,6j and ( |7^ illustrates the 
importance of the bias at moderate p± . 



Figure 3: The prediction for N^g, obtained 
by subtracting from Nq^g the quark contribu- 
tion Nqq, changes significantly if the bias in 
the qq-ferm is neglected. The figure shows 
the effect for y^s^ = 60 GeV, with N^^g 
given by Eg l^^a). 



implies that A^^gg < 2Nf^q for accessible energies. Consequently, the result for Nq^g grows with 
the assumed value of Yq + Y^. 



While the bias is not serious for the gluon term in Eq (|T2D, it is more important for the 
qq term. Focusing on events with comparatively large values of p±, where the bias is less 
essential, and using the assumption in Eq (|l^a), we arrive at the result of ref |jl8[: 



^qqg(^,PiLe) = Kq(Sqq) + ^^ggipl^M^ + (15) 

The bias is formally of order a^, and is here taken into account by the factor (1 + 0{as))- 
The result of this expression, neglecting the 0{as) term, is also shown in Fig ^ The effect 
of the bias corresponds to less than one charged particle for p^cut larger than ~ lOGeV, but 
becomes much more important for smaller p^cut- values. 

An alternative way to express this result is the effect on extracting A'gg from data for A^qqg, 
as illustrated in Fig |. A^gg can be extracted by subtracting the biased quark multiplicity 
A^'qq(Lqq, klu) from iVqqg, hcrc assumcd to be described by Eq (0a). Neglecting the bias in 
the subtracted iVqq term gives a significantly different result. The relative effect of the bias 
is in this case larger, and it exceeds 20% for p± < 15GeV. Furthermore, to get a reliable 
result for Ngg, the relevance of subleading terms in the biased quark multiplicity needs to 
well understood. For the solid line in Fig H, the MLLA relation in Eq ([T0|) is used to subtract 
the qq contribution from the total multiplicity. Instead using the LLA relation in Eq (P) 
would give a prediction for A^gg which is about three charged particles higher for most values 

of P±cut- 

Although the effect of the bias is very important for small p±, we also see from Figs ^ and ^ 
that it can be neglected for large p_L-values, where, thus, the results in ref |0 and Eq ( [T3| ) 
can be safely used. This implies e.g. that the bias is negligible in gluon systems defined 
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as the hemisphere opposite to two quasi-coUinear quark jets, thoroughly investigated by 
OPAL ii. 

It would be very interesting to compare the results in Figs |^ and ^ to experiments. Experi- 
mental data on A'^qqg can be directly compared to the Monte Carlo or MLLA results in Fig ^ 
Data on the difference Ai'qqg — iVqq can be compared either to the predictions in Fig ^ or to 
experimental results for Ngg obtained through one of the methods described in ref We 
have compared the results in Fig |^ with MC simulations, where the p± scale is determined 
by the Durham cluster algorithm. The MC results (not shown) indicate that an analysis 
based on jet reconstruction is accurate enough to illustrate the effects the bias, but perhaps 
not to distinguish between the assumptions in Eq (|T^ ) and We also note that the 

effects described here may have a phenomenological impact on the recent analysis of N^^g ^ , 
which employs the two-scale dependence. 



6 Conclusion 

A series of subtle effects influence an analysis of the difference between quark and gluon jets 
in a real life experiment. In this letter we discuss and clarify effects associated with 

• the definition of p±, 

• the bias from restrictions on subjet pj_, 

• the problem that infrared cutoffs cannot be uniquely defined from perturbative QCD. 

We also demonstrate the impact of these effects on the analysis of three-jet events in e~^e~- 
annihilation. 
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